clc;
clear;
a=0:0.1:4;
f=@(x)exp(-a*x.^2).*sin(a.^2*x);
I=integral(f,0,inf,'RelTol',1e-20,'ArrayValued',true)
% RelTol用来控制精度
% ArrayValued用来确定被积函数是否返回数组值

plot(a,I)